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Abstract 

We study the question whether a possible metastable vacuum state is actually popu- 
lated in a phase transition in the early universe, as is usually assumed in the discussion 
of vacuum stability bounds e.g. for Standard Model parameters. A phenomenological 
(3 + l)-dimensional Langevin equation is solved numerically for a toy model with a po- 
tential motivated by the finite temperature 1-loop effective potential of the Standard 
Model including additional non-renormalizable operators from an effective theory for 
physics beyond the Standard Model and a time dependent temperature. It turns out 
that whether the metastable vacuum is populated depends critically on the value of 
the phenomenological parameter r] for small scalar couplings. For large enough scalar 
couplings and with our specific form of the non-renormalizable operators the system 
(governed by the Langevin equation) always ends up in the metastable minimum. 
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1. The symmetry breaking potential of a field theory at vanishing temperature T may have 
local (energetically disfavoured) minima in addition to the global minimum which defines 
the true vacuum. This leads to metastability if the system sits at T = in a false vacuum 
and the decay of such a false vacuum may occur in two ways. First it is possible to cross the 
classical energy barrier which results from the balance between the released volume energy 
of a single "bubble of the correct vacuum" with the amount of energy stored in the surface 
of that bubble. Vacuum decay sets in after formation of the first bubble of critical size which 
requires however a big enough thermal fluctuation or some source of local energy deposition. 
Alternatively the false vacuum may decay via quantum tunneling and the parameters of the 
theory determine in general which decay mechanism dominates. 

The study of questions concerning the decay of metastable vacua under the influence of 
fluctuations induced by quantum effects goes back to the work of Krive and Linde PJ who 
studied vacuum stability under quantum corrections in the a-model. The issue of vacuum 
stability may also be relevant for the Standard Model (see e.g. the review of Sher || and 
references therein) and for model extensions. Requiring that our electro-weak vacuum is 
the true vacuum leads for a large top mass m t ^80 GeV to lower bounds on the Higgs mass 
0. These bounds can be somewhat relaxed by allowing the electro- weak vacuum to be 
metastable with a lifetime longer than the age of the universe 0, where the dominating 
decay mode is the classical transition via the formation of a critical bubble triggered by 
random collisions of ultra high energy cosmic rays. Due to the experimental lower bounds 
on the Higgs mass metastability is in the Standard Model almost ruled out. It may however 
still be an interesting issue in different model extensions and is certainly also in general an 
interesting question in field theory. 

Thermal effects were included only in a few cases in studies of metastability (see e.g. 0- 
0). Most authors concentrated on transitions at vanishing temperature assuming that 
the universe is already in the false vacuum. This ignores however that the symmetry was 
restored at high temperature in the early universe and we will discuss in this paper the 
question how likely or unlikely it is that the system falls into a metastable vacuum with 
cosmological boundary conditions. Given the physical situation realized in the early universe 
(or more generally in a system that is cooled through a phase transition), we calculate thus 
for a theory exhibiting symmetry breaking with metastable vacua at T = the probability 
to end up in the false vacuum in the first place. 

To address this topic one needs to study the response of a given field theory to a changing 
temperature which is obviously a highly non-trivial taskFJ. A calculation from first principles 
would require to self-consistently identify a heat bath, to calculate the coupling between the 

1 Since the temperature in the situation considered here explicitely depends on time, one can not take 
recourse to well established approaches like Langers theory of bubble nucleation. Langers description 
may however be used in the late stages of the evolution when the time-dependence of the temperature 
is dynamically no longer important and if the potential barrier seperating the vacua is sufficiently high such 
that a saddle point expansion about the top of the barrier is a good approximation S . 
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bath and the system and solve the equations of motion for the macroscopic expectation val- 
ues that describe the system as the temperature of the heat bath is changed. All this would 
have to be done in the framework of non-equilibrium field theory. The conceptual founda- 
tions for such a study are in principle given (at least for a scalar theory) by Morikawa || 
and Gleiser and Ramos flO| (c.f. also Boyanovsky et. al. Jill). Using the CTP-formulation 



of field theory by Schwinger and Keldysh they discuss the (perturbative) formulation 



of an effective equation of motion in a real Ay -theory which turns out to resemble the well 
known Langevin equation in the form 

n<p + V'(<p)+ri(<p)<p = t{<p) , (1) 

albeit with rather complicated dissipative and noise contributions. 

This makes contact with a number of numerical studies that use equations of the type given 
in eq. ([!]) to analyze the decay of false vacua at finite temperature [fL4H-||21]|. In |20| and ||2T 



for example, (3+l)-dimensional theories are studied using an equation of the form 

n<p + U , (<p,T)+r } <p = Z, (2) 

where the used effective potential U(<p, T) is motivated by the high temperature result found 
in the Standard Model written as 

^,T) = ^(r 2 -r 2 V-|r^ + ^ 4 , (3) 

where r\ and £ are related by the fluctuation - dissipation relation 

(^,m?',t'))=2 V T6^\x-x') . (4) 

The authors studied the system at a fixed (namely the critical) temperature and choose r] 
to be 1, arguing that the value of 77 would not affect the final (thermodynamical) state of 
the system, but only the approach to equilibrium. Here we will use a Langevin equation 
of the form @ to model the dynamics of a scalar theory with a symmetry breaking phase 
transition (at some critical temperature T c ) and with metastable vacua for < T( <^T C ) in 
order to calculate the probability of metastable states. As mentioned above the question 
is whether the metastable states decay during cooling by thermal activation long before 
the temperature approaches zero or if the system populates the metastable states with 
reasonable probability which would be the starting of the usual vacuum stability discussion 
in the Standard Model. 

2. We first discuss the choice for the effective potential U((p,T) appearing in equation 
(Q) and for the time dependence of the temperature. Even though we will not aim at 
a discussion of the full Standard Model, embedded in some specific extension at a given 
scale A which renders the physical vacuum metastable, it makes sense to stick as close as 
possible to the Standard Model potential. For reasons to be discussed below, we will only 
be interested in theories with a relatively small scalar self-coupling. We chose therefore as 
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a first contribution to U((p, T) the 1-loop effective potential of the Standard Model in the 
form (c.f. fI3|) 
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where gs,F count the number of degrees of freedom of the corresponding fields and, taking 
into account only the SU(2) gauge bosons and the top quark, 
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with the vacuum expectation value a = 246 GeV. 

The influence of new physics from an extended model at some scale A can be parametrized 
by suitable higher dimensional non-renormalizable operators. For a specific extension these 
operators can be obtained as usual by connecting tree graphs which contain light scalars by 
propagators of heavy degrees of freedom. For illustration one can assume heavy scalars with 
masses M 2 ~ A 2 + T 2 which interact with the light scalars. To obtain the leading contri- 
bution to the six point function one connects two quartic vertices which contain one heavy 
scalar. For small external momenta this results at T = in an effective six point function 
for the light scalar of the form yj 6 /A 2 . Including temperature effects results approximately 
in ip 6 /(A 2 + T 2 ) and for small temperatures we can finally expand in T/A. 

For if <C A the discussion can be restricted to the leading higher dimensional operators. 
While contributions to the two and four point functions can be absorbed into the definition 
of the usual scalar parameters the leading new terms occur in the 6- and 8-point functions. 
This procedure yields thus an additional contribution to the potential of the form 
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which is also heuristically suggestive. For the present study ge and #8 should be viewed 
as couplings which parametrize new physics. The common T-dependent factor has the 
following properties: it vanishes for T — ► A where the theory is not reliable any more and 
it becomes unity for T — > as in the usual zero temperature field theory. 
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Finally, the interesting temperatures are of the order of the weak scale, T ~ (9(100) GeV. 
Taking the time dependence of the temperature to be governed by the (adiabatic) expansion 
of the universe according to the Friedman equations (curvature effects can be neglected 
at these temperatures to a very good approximation) we have for a radiation dominated 
universe with thermodynamical equation of state p = 3p that 

T oc r 1 ' 2 . (8) 



Using Gaussian white noise obeying the fluctuation - dissipation relation (f|) the only free 
parameters are now the scalar self-coupling A (resp. the Higgs mass m H ), the couplings g 6 , 
<?8 and the scale A, as well as the viscosity coefficient rj. We consider different values of A, 
taking A between A = 0.03 (since for the physical Higgs mass one has m 2 H = (2A + 12B)a 2 
with B as in eq. @ which is B ~ —0.0045 if the measured masses are used, A has to be 
larger than 0.028 for m# to be real) and - for reasons to be discussed below - A = 0.044, 
corresponding to m# between 18.4 and 45.1 GeV. For every value of A, we choose values of 
g 6 and g 8 such that the full effective potential U(ip,T) = U(^(ip,T) + U^)(<p,T) for T = 
has in addition to the physical minimum at a & second, deeper one at a' close to the scale 
of new physics A, which we choose to be at A = 4a. This situation is reminiscent of the 
scenarios considered in connection with the stability bounds as discussed above. Finally r\ 
is varied for each A in order to identify possible different dynamical regimes. 

3. Before discussing the results, let us give some details of the numerical procedure used. 
The full details may be found in [2^]. We want to perform numerical simulations in our 
simplified model inspired by the Standard Model in an early universe context. Therefore the 
system is rewritten in dimensionless quantities and we adopt the choice of by writing 
explicitly 



x := o}l 2 T2X 



a := a 3 ^ 4 a 



g e := a 1/2 g 6 



t 



a l ' 2 T,t 



= T := T/T 2 

A(t) := gT 1//2 A(t) 
g 8 := a~ 1/2 g 8 



where a = 2D and 
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Omitting the twiddles, in the dimensionless system eq. (0) has the form 



d 2 X 



with the fluctuation - dissipation relation 

(£(x, t)£(x', f)> = 2rfiS(t - t')5 (3) (x - x' 
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Instead of evaluating the full expression for the potential (|5|) at any temperature we will 
for simplicity rely on the high and low temperature expansion (HT and LT) [[13] in the 
appropriate regions of m/T . It is well known that both the high and the low temperature 
expansion work surprisingly well even for m/T of order 1. Taking into account a field inde- 
pendent contribution to the HT-potential (which is often ignored) HT and LT can therefore 
be matched in a reasonably smooth manner at m = T. Even though there are different 
masses entering the potential we may exploit the fact that all masses are proportional to 
the scalar expectation value and use the top mass to decide upon which of the two regimes 
is appropriate. The potential of eq. @ can now be written as 
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where K u (x) are the modified Bessel functions of the second kind and of the order v. Since 
we need later on the first and second derivatives of U, it is more appropriate to use K 2 
explicitely instead of the more familiar exponential representation. The choices for g§ and 
g 8 depending on A will be given below. Note the explicit time dependence of the potential 
appearing in eq. (JTT|) which is a consequence of the t-dependence of the temperature as 
given by eq. (||). 

This equation is analytically not solvable and we use therefore a discretization scheme for a 
numerical integration. Consider a hyper-cubic lattice in 3+1 dimensions with N x points in 
each spatial direction separated by the spacing 5x and N t points in the temporal direction 
separated by the spacing St. Then Xj >n is the value of X at the lattice point (xi,n) where 
i = (ixjiyjiz)- For sufficiently small spacings the derivatives can be written as differences 
whose exact form depends on the discretization scheme. Like in |2IJ , we use the well known 
staggered leapfrog scheme which is of second order accuracy in time. Next the second order 



differential equation ( IT ) is rewritten as a pair of two first order equations for X and X. 
Introducing additional lattices for X which are shifted by St/2 (i.e. X i ri+1 /2 and Xi^yz) 
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and applying the substitution X itn+ i = X^ n + St • X ijn+1 / 2 yields 

1 
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The noise term becomes in discretized form 



2r]9 



and thus the fluctuation - dissipation theorem now reads 



(14) 



(15) 



(16) 



where the Q i>n are Gaussian distributed random numbers of unit variance. The initial 
conditions are chosen to be 



X„ 



i,0 



and Xj 



'1/2 



o , 



;i7) 



which corresponds to a homogeneous system without fluctuations supposedly describing a 
rapidly and adiabatically cooled (i.e. quenched) state of the early universe before the onset 
of the electro-weak phase transition. We also impose periodic boundary conditions to the 
finite lattice, i.e. compactification of the spatial dimensions. The generation of unwanted 
long range correlations is avoided due to the uncorrelated noise terms. 

The simulation requires a few more words of explanation. Since the temperature evolves 
with time we have to fix the initial and the final temperature, 9 ini and 9^ n , respectively, 
which are connected by 



,1/2 _ a ,V2 



yielding thus the number N t of time steps St, 
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(18) 



(19) 



We choose for all calculations 6> ini = 1.5, well above the onset of the phase transition. 
The final temperature is to a large extent free; one just has to take care that all interesting 
dynamics has happened already when it is reached. 9& n should preferentially not be too small 
since N t depends strongly on the ratio 9^/9^. Usually, 9 &n = 0.15 is a good compromise, 
together with St = 0.1. The temperature 9{t) and thus the potential U(X,9(t)) with its 
first and second derivative are calculated at every time step of the simulation. All extrema 
of U are searched and it is determined whether they are minima or maxima. This fixes the 
number of phases which may be one, two or three in our model. Note that for 9 > the 
second (or third) minimum may be at field values which are larger than the cutoff A. 
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A 


9e 


9s 


m H (GeV) 


o' (GeV) 


^7crit 


0.030 


0.003941 


0.063056 


18.4 


948.2 


0.69 


0.031 


0.003892 


0.062279 


21.4 


941.6 


0.60 


0.032 


0.003578 


0.057254 


24.1 


958.2 


0.54 


0.033 


0.003558 


0.056929 


26.4 


948.1 


0.45 


0.034 


0.003549 


0.056791 


28.7 


938.1 


0.40 


0.035 


0.003285 


0.052557 


30.7 


951.2 


0.34 


0.036 


0.003290 


0.052645 


32.6 


937.4 


0.27 


0.037 


0.003054 


0.048859 


34.4 


949.2 


0.21 


0.038 


0.003066 


0.049053 


36.1 


933.0 


0.15 


0.039 


0.002806 


0.044899 


37.8 


949.6 


0.12 


0.040 


0.002818 


0.045082 


39.3 


930.6 


0.10 


0.041 


0.002620 


0.041916 


40.9 


940.6 


0.08 


0.042 


0.002436 


0.038978 


42.3 


948.8 


0.06 


0.043 


0.002400 


0.038404 


43.7 


933.3 


0.01 


0.044 


0.002227 


0.035640 


45.1 


941.3 





Table 1: Potential parameters (c.f. eq.s (|5|)-(0) and (||)) and corresponding critical values of ?7 cr it 

4. In order to discuss the results of the simulations we need to specify first values for 
the parameters of the potential. As described above, we assume the non-renormalizable 
operators introduced in eq. (^) to be generated by new physics with characteristic scale 
A = 4cr ~ 1 TeV. The corresponding couplings will be chosen such that for any given value 
of A the absolute minimum of the full potential is at a 1 ~ A for vanishing temperature. 
Finally A is varied as mentioned above between 0.03 and 0.044. Table 1 shows the values of 
the couplings A, g 6 and g 8 together with the corresponding values of m H and the position of 
the stable minimum at T = 0, a'. The corresponding potentials at vanishing temperature for 
the parameters of table 1 are plotted in figure 1. Note that height of the barrier between the 
metastable minimum at a = 246 GeV and the stable minimum at a' increases with growing 
A. The system is prepared with X = X = at At = 0, and the temperature is taken to be 
1.5 x T 2 initially. The temperature is then taken to depend on time according to eq. (§) and 
the Langevin equation is solved as described above. A typical result of a solution of equation 
(Q) with the above assumptions on the potential and for different values of r\ is displayed in 
figure 2, where we plot the first and second moment of the dimensionless field, namely the 
mean value of the field (X) and the second moment AX = J ((X — (X)) 2 ) as functions of 
t for A = 0.03 and different values of rj. After an initial phase where the behaviour of the 
system is rather independent of r\ and where it evolves essentially homogeneously (AX is 
very small), the behaviour of the first moment shows a strong dependence on the damping 
parameter 77: For small 77 the system reaches the absolute minimum relatively fast and in 
a more or less homogeneous fashion, performing damped oscillations around X = a'. As 
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X=0.030 



Figure 1: The dimensionless potentials U as functions of the dimensionless field X for vanishing 
temperature with #6,8 as given in table 1. With increasing values of A the height of the barrier 
between the metastable and the true vacuum gets larger. 

7] gets larger, the system remains close to the metastable minimum at a m a' /4 for longer 
times and the transition to the absolute minimum, though still occuring for 77 < 0.7, takes 
longer and longer (note how the second moment deviates from for longer times - the 
system is highly inhomogeneous during the transition). For r\ = 0.7, the behaviour changes 
dramatically: The system is no longer able to cross the potential barrier and to populate 
the true vacuum, but remains in the metastable minimum with only very small fluctuations. 
The region of 77 around the critical value r] crit = 0.69 is resolved in the lower panel of figure 
2. We perform the same procedure for the sets of parameters as given in the table, always 
identifying the critical value of r\ as long as it exists. From table 1, one notices that 77 crit 
is a decreasing function of A. In figure 3 we display the first moment for a relatively large 
coupling A = 0.42. The critical 77 turns out to be 0.06 in this case. The dynamics of the 
system differs from the case displayed in figure 2 mostly by more pronounced oscillations 
around the false vacuum state for small At. These may be understood by noting that first 
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Figure 2: The first and second moments (X) (upper curves) and AX (lower curves) for A = 
0.03. In the lower panel the dependence on rj around r) cr \t = 0.69 is resolved (note the different 
timescales) . 
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Figure 3: The first moment for A = 0.042. 



the values of r] displayed are small compared to the ones shown in fig. 2, i.e. the system is less 
damped, and that second the curvature of the potential at the metastable minimum is larger 
for larger A (c.f. figure 1). By a simple kinematic analogon one expects more pronounced 
oscillations. Finally, for A = 0.044, the critical value of rj vanishes and the system is no 
longer able to reach the true vacuum by thermal activation, irrespective of the value of rj. 

We have checked finite size effects and found them to be small for the lattice size chosen. 
In particular the values for rj CT - lt found on a lattice with 64 3 spatial sites are identical up to 
±0.01 with the ones found on a 50 3 lattice and given in the table. We have also performed 
calculations with time steps smaller than 0.1 as used for the results given above. Thus 
e.g. the results for the moments at A = 0.03 and rj around 0.69, where one expects the 
errors of discretization to be most pronounced, are identical within the numerical precision 
for 5t = 0.1, 0.05, and 0.03. Too large values for St result of course in large deviations 
from the behaviour found above. We thus believe that the numerical error on i] crit (for the 
potential used) is ±0.01 and thus well under control. 

5. To conclude, we have discussed the dynamics of the population and decay of metastable 
vacua in an environment with a time dependent temperature. We have done this using a 
phenomenological Langevin equation for a scalar field with an effective potential motivated 
by the Standard Model embedded in some model with characteristic scale A ~ 1 TeV, 
leading for vanishing temperature to the existence of a stable vacuum state at <p — a' ~ A. 
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This is the situation assumed for the derivation of bounds on Standard Model parameters 
from vacuum stability as considered in the literature . We have discussed the dependence 
of the dynamics on the (phenomenological) parameter 77, governing dissipative effects. In 
particular to answer to the question whether the theory, starting at high temperatures in 
the symmetric phase, populates the (now metastable) physical vacuum depends critically on 
this parameter. In the framework studied here, using an equation of motion of the form @ 
with Gaussian white noise £ related to 77 by (||) and a potential given by we find that 

the metastable vacuum is indeed populated for any value of 77 if A ^0.043, corresponding to 
mil ~ 43.7 GeV. We have shown that for smaller values of A there exists a critical value of 
the viscosity parameter, 77^, such that for 77 < 77^ the system reaches the stable vacuum 
by thermal activation before T = 0, whereas for 77 > r) cr i t the system will populate the 
metastable vacuum (being closer to tp = in our case) and is not able to cross the potential 
barrier by thermal activation. 

Clearly the results collected in table 1 depend both on the specific form of the Langevin 
equation and on the effective potential. However, they highlight the fact that one may under 
certain conditions hope to learn something about the dynamics of phase separation without 
detailed knowledge of the parameters of the effective equation of motion. Thus, while for 
small A the late time behaviour of the model studied here crucially depends on whether 
77 is larger or smaller than ?7 crit , for large enough quartic coupling the value of 77 is indeed 
irrelevant. For a given theory it might thus suffice to calculate the effective potential and its 
dependence on the temperature if the parameters of the theory are in a range where 77 C rit = 0. 
With the specific assumptions made in this work our results thus justify the usual discussion 
of metastability in the Standard Model. On the other hand we also note that for potentials 
with multiple minima the final state of the system - neglecting quantum effects, i.e. in a 
treatment on the basis of Langevin type equations - is not generally independent of 77 if the 
system is "quenched", i.e. the temperature is time dependent. This fact should motivate 



more efforts along the lines of M, [10] and [11] in order to improve our understanding of the 



derivation of phenomenological equations of motion and their parameters from field theory. 
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